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Abstract 

In this paper we simultaneously explain the excesses of the 750 GeV diphoton, muon g-2 and h — 
LT in an extension of the two-Higgs-doublet model (2HDM) with additional vector-like fermions 
and a CP-odd scalar singlet (P) which is identified as the 750 GeV resonance. This 750 GeV 
resonance has a mixing with the CP-odd scalar (A) in 2HDM, which leads to a coupling between P 
and the SM particles as well as a coupling between A and the vector-like fermions. Such a mixing 
and couplings are strongly constrained by T — py, muon g-2 and the 750 GeV diphoton data. We 
scan over the parameter space and find that such an extension can simultaneously account for the 


observed excesses of 750 GeV diphoton, muon g-2 and h — pt. The 750 GeV resonance decays 


in exotic modes, such as P > hA, P — HZ, P ~ HA and P —^ W-HT, and its width can be 


dozens of GeV and is sensitive to the mixing angle. 


PACS numbers: 12.60.Fr, 14.80.Ec, 14.80.Bn 


I. INTRODUCTION 


Very recently, the ATLAS and CMS collaborations have both reported an excess of 750 
GeV diphoton resonance [1], with a local significance of 3.60 and 2.60 respectively. Combin- 
ing the 8 and 13 TeV data, the production cross section times the branching ratio is around 
4.47+1.86 fb for CMS and 10.642.9 fb for ATLAS [2]. However, there are no excesses for 
dijet [3], tt [4], diboson or dilepton channels, which gives a challenge to possible new physics 
explanations of the 750 GeV diphoton resonance [2, 5-13]. 

In addition, the CMS has reported a 2.40 excess in the lepton-flavor-violating (LFV) 
Higgs decay À — pr (here h is the 125 GeV SM-like Higgs), i.e., Br(h — uT) = (0.84*032)96 
[14], while the ATLAS data is Br(h — ur) = (0.740.62)% [15]. This excess can be explained 
in the general two-Higgs-doublet model (2HDM) with LFV Higgs interactions. Also such 
a model can give a sizable positive contribution to the muon anomalous magnetic moment 
(muon g-2) and accommodate the long-standing anomaly [16-18]. 

Attempting to simultaneously explain the excesses of 750 GeV diphoton, h — ur and 
muon g-2, we in this work introduce additional vector-like fermions and a CP-odd scalar 
singlet (P) to the general 2HDM. The singlet P is identified as the 750 GeV resonance, 
which has a mixing with the CP-odd scalar (A) in the original 2HDM. Therefore, the model 
can lead to the P couplings to SM particles and the A couplings to vector-like fermions. 
In addition to the 125 GeV Higgs and 750 GeV resonance data, the LFV Higgs decay 
T — uy can give strong constraints on the couplings and mixing. The dominant decays of 


the 750 GeV resonance can be some exotic modes, such as P > hA, P — HA, P > HZ 


and P — W-H'*. Considering various relevant experimental constraints, we examine the 
diphoton production and decay of the 750 GeV resonance, as well as muon g-2 and h > pr. 

Our work is organized as follows. In Sec. II we introduce additional vector-like fermions 
and a CP-odd scalar singlet to the 2HDM. In Sec. III we perform numerical calculations 
and discuss the muon g-2, h > uT and the diphoton production and decay of the 750 GeV 


resonance in the allowed parameter space. Finally, we give our conclusion in Sec. IV. 


II. MODEL 


We introduce a CP-odd scalar singlet field P to the general 2HDM with the assumption 
that P does not develop a vacuum expectation value (VEV). The Higgs potential is given 


by [18, 19] ; 
£M 


1 Pt — ip Polo, + h.c., (1) 


1 
V= VonpM + 5h PS + 


with 

VonpM = (bib) + io (9355) + used; + he. 
X (910;)? + Ay (0105)? + A(DiD1)( Di Do) + ADi) (51) 
- (812)? Pi h.c.| $ [As(®{1)(®} 2) + he. 


E I(9195)(9[:) + hc. s (2) 


In the Higgs basis, the ®; field has a VEV v =246 GeV, and the VEV of ®, field is zero. 


The two complex scalar doublets with hypercharge Y = 1 can be expressed as 
Gt H+ 
(OR = i , " o. = 1 l . (3) 

73 (v + pi + iGo) 73 (P2 + iAo) 


The Nambu-Goldstone bosons G? and G* are eaten by the gauge bosons. The physical 


CP-even Higgs bosons h and H are the linear combinations of pı and po: 


Pi cosa sina h 


p2 —sina cosa A 


where tan 2a = 2Agv?/(m?5, — M?) with 


1 
Mary = 2d10?, Moo = mays zig v (334 + As). (5) 


'The masses of two CP-even Higgs bosons are given as 


1 
mpg = 5 om + mios F V ant — Moog)? + 4204 | . (6) 


The field HT is the mass eigenstate of the charged Higgs boson, and the CP-odd Higgs field 


Ao has a mixing with P: 


Ap cos? —sind A 
Py sin 0 cos 0 P 


where tan 20 = 2uv/(mA, — m?) with 


1 
mA, = mus + (M — As). (8) 


The masses of two CP-odd scalars are given as 


1 
MAP = E má, + mph F "TS — mà)? + 4u?v?|. (9) 


The 750 GeV Higgs boson P couplings to other Higgs bosons and gauge bosons as 


1 
PAh: cesev [(As + À4 — 2A5)Ca — À7Sa] — a; (mA — m2)s4gCa, 


1 
PAH : cesev [(A3 + A4 — 2À5)8a + Ave] — — (mà, — m2) 84680, 


Av 
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The general Yukawa interactions of the SM fermions are given by 


-L = Yu Qr Pi ur + yaQ, Pi dr + y Ly Per 
+Q, Pour + DO: o. dg + "am s €n + h.c., (11) 


where QT = (ur ,dr), LT = (vz, Iz), $,» = ir >, and Yu, Ya, Ye, p", p! and p* are 3 x 3 
matrices in family space. 

Also, we introduce a singlet quark with 2 electric charge and multiple singlet leptons. 
The Yukawa interactions of vector-like fermions are written as 


-L=mr T T +iyr Py T 5s T - V (mii Li Li t y Po Li ys Li). (12) 
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Then we obtain the Yukawa couplings of the neutral Higgs bosons: 


ij m aC ij V2 a) ij Ü aC ij V2 CE 
jl, " 
jam ci a (for u), Yai = ign (for d, £), 
f. i 


YPij = iis (for u), ypj = —i—4 59 (for d, 2), 


YATT = lyrSe, UAL,L, = İYL; 86; 


YPTT = lyrCe, UPL;L; = lyr;Co- (13) 


For the diagonal matrix elements of o", p^ and pf, we take 


u v2m? des V2mi k £ V2m! n 
v 


Pii " Ku, Pi d: Pü = $ 


t (14) 


which corresponds to the aligned 2HDM [20]. We assume that fis and pis are nonzero, and 
other nondiagonal matrix elements of o", p^ and pf are zero. 

The vector-like quark is introduced to make the 750 GeV Higgs singlet to be produced via 
the gluon-gluon fusion process. However, the vector-like quark can also enhance the cross 
section of gg — A, which will be constrained by the experimental data from the ATLAS 
and CMS searches. Therefore, we expect that the vector-like leptons play the main role in 
enhancing the 750 GeV diphoton production rate. The decay P — yy can be enhanced by 
the vector-like leptons, and its amplitude is proportional to the couplings and the square of 
electric charge. Here we do not discuss the electric charge and coupling of every vector-like 
lepton as well as the quantity of vector-like leptons in detail; instead we focus on the total 


contribution of vector-like leptons, which depends on 
Yz = Da (15) 
where L; denotes the i-th vector-like lepton. 
III. NUMERICAL CALCULATIONS AND DISCUSSIONS 
A. Numerical calculations 


In our calculations, we scan over the parameters in the following range 


—0.06 < Sa < 0.06, — 0.3 < sọ < 0.3, 

0.05 < pur = Pry < 1, — 50 < ke < 50, 

0 < Yz < 50, 0 < À3, À7 < 47, 

200 GeV < my < 450 GeV, (16) 


and fix 


mp = 125.5 GeV mp — 750 GeV, mye =m, = 500 GeV, 


mr = 400 GeV, mr = 700 GeV, yr=2.0, 
Ku = Kq = O. (17) 


During the scan, we consider the following experimental constraints and observables: 


(1) Precision electroweak data. According to the expressions for the oblique parameters 
S, T and U in the 2HDM [21], for —0.06 < s, < 0.06 and c, 1, the expressions in 


this model are approximately given as 


5 = mm, [Ac Fs (mz, mir, m) T c, sgFs(mz, mir, mp) a Fs(mZ, mis, mx) ’ 
1 : 
T = 16 2 2 cc F'r(miy, mA) = ci saFr(my, mp) T c2 Fr(mis, mj) 
TMy Sy 


+ e Fr(mijs, mA) + sS Fr(miys, mp)] , 


U = 2 [A Fs(my, mys, Mma) — 2Fs(myy, mi Tg ) 
T TI 


+G Fs(myy, Mass m) + SES (my, Ms, my] 


1 
ns [sco Fs(mz, mir, mA) + css Fs(mz, mg, mp) 
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— Fs(mZ, mj Sm )] ; (18) 
where 
1 ab a 
Fr(a, b) 2s 5 (a b) = gh log(7) Fs(a, b, c) ms B»»(a, b, c) P3 B22(0, b, c) (19) 
with 
Bodas ites [ex (X — de) 
22(a,b,c) = 7 |b+c— 3a 2 Jj, x X log ie), 
X = bx + c(1 — x) — ax(1 — x). (20) 
Here we require [22] 
S = —0.03 + 0.1, T= 0.01 40.12, U = 0.05 0.1 (21) 


(2) The 125 GeV Higgs data. For —0.06 < Sa < 0.06, ku = Ka = 0 and —50 < ky < 50, 
the 125 GeV Higgs couplings to the gauge bosons, up-type quark and down-type quark 
are very close to the SM values, but the coupling to 7T can have a sizable deviation 
from the SM value. The signal strength of 77 channel is à,, = 1.41793, from ATLAS 
[23] and fi. = 0.89*031 from CMS [24]. We require 0.33 < fi, < 2.21 and such a 
bound will give strong constrains on s, and ke for which the absolute value of the 


coupling of the 125 GeV Higgs and 77 is around the SM value. 


(3) 


Non-observation of additional Higgs bosons. For —0.06 < sa < 0.06 and Kku = Ka = 0, 


the cross sections of H and HT at the collider are very small, and hence H and HF can 
be hardly constrained by the current experimental data from the ATLAS and CMS 
searches. The pseudoscalar A can be produced via the gluon-gluon fusion process with 
vector-like quark loop, and the decay A — yy, A — yZ and A — ZZ can be enhanced 
by the vector-like quark and leptons at one-loop level. For m4 =500 GeV, we impose 


the following relevant bounds at the 8 TeV LHC [25-29] 


FR < 6 fb, Rzg < 45 fb, Rzy « 6.8 fb, Raz < 60 fb, Re, < 26 fb. (22) 


The 750 GeV resonance data. The 750 GeV Higgs singlet P can be produced via 
the gluon-gluon fusion process with vector-like quark loop, and the decays A — yy, 
A — JZ and À — ZZ can be enhanced by the vector-like quark and leptons at one- 
loop level. Due to the mixing with A, the 750 GeV singlet can decay into the 5M 
particles, such as hZ and 77. For the 750 GeV Higgs singlet, we impose the following 
bounds at the 8 TeV LHC [25-29] 


Ry < 2 fb, Rzz < 12 fb, Rzy <4 fb, Raz «19 fb, Rz < 12 fb. (23) 
At the 13 TeV LHC, we require the 750 GeV diphoton production rate as 
2 fb < Ry < 10 fb. (24) 


The data of Br(h — ur). The branching ratio of h — ur is given by 


(Par + Pou) ma 
167T y, 


where I; is the total width of 125 GeV Higgs. To explain the h — uT excess reported 


Br(h > pr) = (25) 


by CMS within 20 range, we require 


0.1% « Br(h — pr) < 1.62%. (26) 


The muon g-2 data. The dominant contributions to the muon g-2 are from the one-loop 


diagrams with the Higgs LFV coupling [30], 


5 MMs Pur Pry sa (log 73 = 2) c; (log m2 
"boss 167? m2 ' m 
h H 


m2 m2 
— T — z f 27 
= = (27) 


The muon g-2 can be also corrected by the two-loop Barr-Zee diagrams with the 
fermions loops, W and Goldstone loops. Using the well-known classical formulates 


[31], the main contributions of two-loop Barr-Zee diagrams in this model are given as 


amy, 5 
dau = E > N? QF Youn Yost Fo(v5o) 
Í &—h H,A,P;f—tbr T bs 


am 23 
+ gas D Youn goww E (two) + Fa (awe) 
o=h,H 
3 m3 
rid (owe) + Im? {Fn (two) — FA(xwo)] | (28) 
W 


where tg = m?5/m2, twe = miy/m2, Qguww = Sa, g9aww = Ca and 


Fo = Fuly) = E f dem egt) for d=h, H) (9) 


2 a(l—z)-y 
F9) = Fa) =i f de — quoa P) (80 
ENS ee Lo rt) 
fer] 0m 


The experimental value of muon g-2 excess is [32] 


õa, = (26.2 + 8.5) x 107". (32) 


The data of Br(r — py). The LFV coupling of the Higgs boson gives the dominant 
contributions to the decay 7 — uy. The branching ratio of 7 > uy is given by 


BR(T — py) = 481?a (|Airo + Aire + Aor]? + |Aino + Aire + Ar|?) 


33 
BR(r > puv,v;) G2, ; (33) 


where Aizo, Aire, Aigo and Are are from the one-loop diagrams with the Higgs bosons 


and tau lepton [17], and 


* 2 
Yorn |,» má  3Y, Vere 
Aizo = ` 2,42 lu TT (ioe 2 = | D (34) 
PONE P Oe 167 mj Wie? 2 6 
ra ud 
A = — 35 
1Lc 19272m2,_ ; ( ) 
Airo = Auro (u$ ru 7? Your, Uó rr © Vo - , (36) 
Aire = 0. (37) 


Here Ag; and A»g are from the two-loop Barr-Zee diagrams with the third-generation 


fermions loops, vector-like fermions loops and W loops [17]: 


Aor = — ` Or [Qy (Re(yo rr) Fu (ye) — ibm(yo rp) FA (@ Fo) 


8r? m,m 
$=h,H,A,P; f=t,b,T,T,Li TCU 
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2 
1— 4s% 2 1- ty \ a 
TT re Ë 
852, l ( wt rw; H(twé wz) 
o 1-t,\ = 3 
+| 7— 3tw — py» FA(xwe,Xwz)- 5 {Fa(twe) + G(zwe)}?|, (38) 
Aor = Aor (u$ TH EE Uo ur; i —> —i) ; (39) 
where T3; denotes the isospin of the fermion, tiy = siy/ciy, v;jg = m4/m3 and 


rwz —mi,/ml, and 


= _ &F p(y) — yFn(x) 


Fa(x,y) = 2 =Y i 2 
dade 6) E M (41) 


The terms in the first two lines of Eq. (38) come from the effective @yy vertex and 
oZy vertex induced by the third-generation fermion loop and vector-like fermion loop. 
Other terms are from the effective oyy vertex and Zy vertex induced by the W-boson 
loop. The current upper bound of Br(r — py) is [33, 34] 


Br(r — py) < 44 x 107$. (42) 


B. Results and discussions 


In Fig. 1, we project the surviving samples on the planes of Ke versus Sa, Pur versus 
mg, Yr, versus Ke, and pur versus Sy. The upper-left panel shows that there is a strong 
correlation between s, and kK, due to the experimental constraints of the 77 channel data of 


125 GeV Higgs. The surviving samples have two different 125 GeV Higgs couplings to 77, 


00 EU 
m, (GeV) 


07? = 


Lol i| 1 
0.04 


FIG. 1: Under the constrains of the oblique parameters and the LHC Higgs data, the surviving 
samples projected on the planes of Ke versus Sa, Pur versus My , Yr, versus Ke and pur versus Sq. 
The circles (green) are allowed by the muon g-2, the pluses (red) allowed by the muon g-2 and 


Br(r — uy), and the bullets (black) allowed by the muon g-2, Br(r > uy) and Br(h > ur). 


and their absolute values are around the SM value. One is the SM-like Higgs coupling with 
the same sign as the coupling of the gauge boson, and the other is the Yukawa coupling with 
the opposite sign to the coupling of the gauge boson for a relative large &;. 

From the upper-right panel of Fig. 1, we see that the muon g-2 favors p,- to increase 
with my. As shown in Eq. (27), the muon g-2 can obtain positive contributions from the 
H loop and negative contributions from A and P loops for p, = p,,. With the decreasing 


of the mass splitting of H and A, the cancelation between the contributions of H and A 


10 


25 30 35 40 


T, (GeV) 
FIG. 2: Under the constrains of the oblique parameters of electroweak data, the LHC Higgs data, 
muon g-2, BR(r > jy) and Br(h — ur), the surviving samples projected on the planes of T p 
versus Sq, Yr, versus I'p and Yr versus sg. Here 2 fb < Ry, < 4 fb for the circles (green), 4 fb 
< Ry < 6 fb for the pluses (red), and 6 fb < R,, < 10 fb for the bullets (black), with R,, denoting 
the 750 GeV Higgs production rate at the 13 TeV LHC. 


loops becomes sizable so that a large p, is required to enhance the muon g-2. From the 
lower-left panel, we see that the upper bound of 7 — py favors a large absolute value of «g 
for a large Yz. The vector-like leptons with a large Yz can sizably enhance Br(r — uy) via 
the two-loop Barr-Zee diagrams, and such contributions can be partially canceled by the 
one-loop diagram for a properly large &; and a proper sign of sg. The lower-right panel shows 
that the experimental data of Br(h — uT) requires p,, to increase with the decreasing of 
the absolute value of sa, and -0.05 < Sa < 0.05 and 0.05 < pur < 0.7 are favored by the 
Br(h > ur), Br(r > py), muon g-2 and the other experimental constraints as mentioned 
above. 

In Fig. 2, we project the surviving samples on the planes of the total width of 750 GeV 
singlet (T p) versus sọ, Yr, versus I'p and YT, versus 59. From the left panel, we can find that 
I'p is very sensitive to sg since the P couplings to SM particles are relevant to sg. The Ip 
value increases with the absolute value of sọ, and reaches 35 GeV for |sg| =0.3. No matter 
how large the Ip value is, the 750 GeV diphoton production rate R,, can vary from 2 fb 
to 10 fb. The middle panel shows that, with the increasing of the total width, Y; becomes 
large enough to enhance Br(P — yy), and further make R,, to be in the range of 2 fb and 
10 fb. For lp = 35 GeV, R > 4 fb requires Yz to be larger than 30. The right panel shows 


that a large absolute value of sg favors a large Yr since sg with a large absolute value will 
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FIG. 3: Under the constrains of the oblique parameters of electroweak data, the LHC Higgs data, 
muon g-2, Br(r — py) and Br(h — pr), the surviving samples with 20 GeV < lp < 40 GeV 
projected on the planes of the widths of the main decay modes of the 750 GeV Higgs versus Às, mg 
and àz. Here T(P — hA) for the circles (green), T(P — HZ) for the triangles (blue), r(P — HA) 


for the pluses (red) and T(P —^ WHF) for the bullets (black). 


enhance the total width of 750 GeV Higgs sizably. 
In Fig. 3, we project the surviving samples with 20 GeV< TP? «40 GeV on the planes 


of the widths of the main decay modes of the 750 GeV Higgs singlet versus A3, my and Az. 
This figure shows that P — Ah, P > HZ, P — AH and P — W-H* are the main decay 


modes, and the decay P — hZ is insignificant due to the suppression of sg. The decay 
P — Ah is sensitive to A3 and increases with A3. The width of P — Ah can reach 20 GeV 
and dominate over other decay modes for A4 = 47. The decay P — HZ is sensitive to my 
and decreases with the increasing of my. The width of P — HZ can reach 10 GeV and be 
larger than those of P — AH and P — WHF for my = 200 GeV. The decay P > AH 
increases with A; and can be larger than the width of P — W-^H* for A; = 10. The width 
of P — W*HF is in between 2 GeV and 4 GeV for my+ = 500 GeV, and not sensitive 


to mg, A3 or Az. With the increasing of m4 and my+, the decay P > hA, P — HA and 
P — W-H* will be kinematically forbidden, which will reduce the width of 750 GeV Higgs 


sizably. 
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IV. CONCLUSION 


To simultaneously accommodate the excesses of the 750 GeV diphoton, muon g-2 and 
h — uT, we proposed an extension of 2HDM with vector-like fermions and a CP-odd scalar 
singlet P, which is identified as the 750 GeV resonance. There is a mixing between the 
750 GeV Higgs and the CP-odd scalar A, which leads to the P coupling to SM particles 
and A coupling to vector-like fermions. In the 2HDM the Higgs bosons have tree-level 
LFV interactions with u — 7, which can be responsible for the excess of h — uT and also 
give sizable contributions to the muon g-2. The 750 GeV Higgs can decay into P — Ah, 
P — HZ, P — AH and P > W-H'*, and its total width is sensitive to sg and can reach 35 
GeV for |s9| = 0.3. Since the 750 GeV Higgs has a large width, the vector-like leptons are 


required to enhance Br(P — yy) to obtain R,, > 2 fb. Meanwhile, such vector-like leptons 
will give sizable contributions to Br(r — uy) due to the mixing of P and A. Therefore, 
the Higgs couplings to 77 are required to be properly large to cancel the contributions of 
vector-like leptons to Br(r — py). Considering the current constraints of the LHC data, 
precision electroweak data and Br(r — py), we scanned over the parameter space and found 
that such an extension can simultaneously explain the excesses of the 750 GeV diphoton, 


muon g-2 and h > pr. 
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